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Abstract. Classical single-mode fiber interferometers, using one fiber per aperture, have very limited imaging 
capabilities and small field of view. Observations of extended sources (resolved by one aperture) cannot be fully 
corrected for wavefront aberrations: accurate measurements of object visibilities are then made very difficult 
from ground-based fiber interferometers. These limitations are very severe for the new generation of interferom- 
eters, which make use of large telescopes equipped with adaptive optics, but can be overcome by using several 
fibers per aperture. This technique improves the wide field imaging capabilities of both ground-based and space 
interferometers. 
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1. Introduction 

In the visible and near-infrared, single-mode optical fibers 
allow efficient transport of coherent light over large dis- 
tances: the transmission losses inside the fiber are typi- 
cally IdB/km in the near-infrared. The throughput of a 
fiber interferometer can be very high, since the number 
of optical elements before the fiber is small, and there is 
no need for numerous beam steering mirrors. Moreover, 
the beam combination can be done using fiber optics cou- 
plers, and does not require high numbers of glass to air 
interfaces. A second interesting property of single mode 
optical fibers is the spatial filtering of the incoming wave- 
fronts. Spatial filtering with single- mode optical fibers 
is superior to pinhole filtering (Keen et al. 2001) and is 
used in interferometry to restore coherence of the tele- 
scope's beams. Single-mode fibers only allow the propa- 
gation of a fundamental mode (LPqi)- The loss of light 
coherence of the entrance pupil due to wavefront aber- 
r ations is then traded for a loss of coupling efficiency 
( ICoudc du Foresto et al. 2000 ). The object visibility mea- 
surements can be corrected for this loss of coupling effi- 
ciency by recording in real time the amount of light cou- 
pled in each fiber (real-time self calibration of the inter- 
ferometer), which greatly improves the accuracy of the 
object visibility measurement. Thanks to these two prop- 
erties, fiber optics offer a reliable and simple solution to 
transport and combine coherent beams in interferometers 
with few apertures, and the recent progress in integrated 
optics will soon allow complex beam combiners to be 
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built for larger interferometer arrays ( Malbet et al 1999 , 
Bergcr ct al. 19991 |Bcrgcr ct al. 200l| ). 

Unfortunately, the restoration of coherence by spatial 
filtering comes at the expense of a very small field of 
view (FOV), which is the size of the diffraction spot of 
the largest telescope in the array. This is a serious lim- 
itation for the new generation of interferometers (Keck, 
VLTI, OHANA) which make use of large AO-equipped 
telescopes: the FOV is very small for such large telescopes. 
Moreover, these interferometers make it possible to ob- 
serve fainter non-stellar sources such as AGNs and YSOs, 
for which the required FOV is often larger than for stel- 
lar observations. In this work, it is shown in §2 and §3 
that spatial filtering reduces the field of view and that 
wavefront aberrations corrupt the object visibility mea- 
surements. This later effect is quantified and discussed in 
§4 and §5. The use of multiple fibers per telescope offers 
an attractive solution to these limitations (§6) and greatly 
improves the wide field imaging capabilities of space (§7) 
and ground-based (§8) interferometers. 



2. Coupling efFiciency and field of view 

2.1. Coupling efficiency 

Coupling of starlight into a single-mode fiber is achieved 
by positioning the fiber head in the focal plane and cen- 
tering it on the PSF of the star. The coupling efficiency, 
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which is the fraction of the starlight coupled into the fiber, 
is given by 



Coupling efficiency of a single-mode fiber 



\IE. 
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J\Ef,\'dAxJ\E, 



011^ dA' 



where Efp is the incident electric field at the head of 
the fiber, in the focal plane, and Eqi is the electric 
field of the fundamental (LPqi) mode of propagation 
in the fiber. Efp is obtained by Fourier transform of 
the electric field in the entrance pupil, and Eqi is very 
well approximated by a Gaussian of constant phase 
across the width of the fiber. The coupling efficiency 
between the LPpi mode an d the Airy fu nction is 0.82 
( ^haklan fc Roddicr 198^ , [Ruilicr 1998] ). Therefore, 
under optimal conditions, no more than 82% of the 
incoming starlight can be coupled into the single-mode 
fiber. Fresnel refiection by the surface of the fiber head 
reduces this coupling efficiency to 78 % unless an anti- 
reflective coating is used on the fiber head. 
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Fig. 1. The coupling efficiency of the light from a point 
source to a single-mode fiber at the telescope focus is a 
function of the position of the source on the sky. 



2.2. Field of view of the fiber 

The coupling efficiency decreases rapidly with the dis- 
tance of the point source to the optical axis of the tele- 
scope. Using equation (1), Shaklan and Roddier (1988) 
have shown that the coupling efficiency, in the case of a 
circular aperture and a fundamental mode approximated 
by a Gaussian, is 

p{a,d,f,uj,X) = -^e V - J 
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Where / is the telescope focal length, d is the telescope 
diameter, A is the wavelength, uj is the 1/e width of the 
fundamental mode, a is the angular position of the point 
source relative to the optical axis, Ji is the first order 
Bessel function and /o is the zero order modified Bessel 
function. We note that this expression is only a function 
of Xi \a\/{X/d) and X2 = (/A)/(wd) : 



p(Xi,X2) = 8 



_o^-2{XiX2f 



J e~^" Io{2X^X2r).h{f-)dr 



(3) 



where r is now dimensionless. p(0,X2) is maximum for 
X2 = 1.401, and this value will be used in this work. We 
can then express p as a function of Xi only, which leads 
to : 



p{a) = 8 e 



-3.923 
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The relation between the amplitude E of the electric 
field coupled in the fiber and the luminosity I of a point 
source which position is ot is then 



E{c.)^^^-^^). 



(5) 



Figure 1 shows the coupling efficiency as a function of | a | , 
in units of A/d, for an unobstructed circular pupil and 
a circular pupil with an central obscuration of diameter 
0.4 X d. In each case, / has been chosen to maximize the 
couphng efficiency for \a\ = 0. At |q;| — X/d, the couphng 
efficiency is 10 times smaller than it is on-axis for a non- 
obstructed circular aperture, and p{0.58X/d) — 0.5p(0). 

Therefore, the field of view of an interferometer using 
one single-mode fiber per aperture to transport the beam 
is limited to X/d. This very limited field of view (FOV) is 
a serious constraint for the observation of many objects, 
such as active galactic nuclei (AGNs) and young stellar 
objects (YSOs). Aperture synthesis using either the Earth 
rotation for ground-based interferometers or a rotation of 
the array for space interferometers could allow reconstruc- 
tion of images with a much larger FOV than permitted by 
the use of single-mode fibers. 

A point of concern is the use of single-mode fibers on an 
interferometer whose telescopes have different diameters 
and/or pupil shapes. The coupling efficiency, as projected 
on the sky, is then different for each telescope, and this ef- 
fect has to be taken into account for image reconstruction 
of objects larger than X/dmax {dmax is the diameter of the 
largest aperture in the array). Residual wavefront aberra- 
tions also produce a similar discrepancy between the fiber 
coupling efficiencies. These effects are studied in the next 
2 sections. 

3. Observation of extended sources with a fiber 
interferometer 

3.1. Instrumental response to a point source 

A general description of fiber interferometers is given in 
Coude du Foresto et al. (1997), and we only recall here 
some of the basic concepts for a 2 telescope interferometer 
(Fig. 2). 
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Telescope i 




with the baseline vector A; = fci — fc2. If a 50/50 fiber 
splitter is used, there are in fact two interferometric out- 
puts (Fig. 2) : 



hit) 



and 



El 
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El 
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E1E2 cos{a ■ k + ij{t)) 



E1E2 cos{a ■ k + ij{t)). 



(10) 



(11) 



Telescope 2 



"1 2 
Interferometric outputs 



The OPD between the two beams is usually modulated 
{il^{t) is not constant), and the temporal modulation of / 
obtained by this changing OPD leads to the measurement 
of El + El (mean level), E1E2 (fringe amplitude) and 
(j) = OL-k (fringe phase). For example, most interferometers 
introduce a time-dependent OPD variation to scan the 
fringe packet and extract the fringe parameters in each 
scan. 

The phase of the interferometric signal, (j) — a ■ k, can 
only be referenced through the simultaneous observation 
of the object and a reference star or by using phase closure 
techniques, which requires at least 3 apertures. 



Fig. 2. Optical la yout of a 2 telescopes fiber int erferome- 3 2. Instrumental response to an extended source 



ter. Adapted from Coude du Foresto et al. 1997, The light 
of each telescope is coupled into a fiber in the focal plane. 
Before interferometric combination in fiber coupler X, the 
coupling efficiency in each arm of the interferometer is 
measured (photometric outputs Pi and P2) thanks to the 
fiber couplers Yi and 12- 



A fiber interferometer cophased (zero optical path- 
length difference (OPD) from the source to the detector) 
for a direction Tq on the sky couples the electric field of 
the observed source on the telescopes' focal planes into 
the fibers. With a 2 telescope interferometer observing an 
unresolved point source at a direction r, the electric fields 
El and E2 in the fibers, as a function of time t, are given 

by 



The expression of the interferometric signal for an ex- 
tended source, of light distribution I(q:), is obtained by 
integration of equation (9) over a : 

I{t) = / El{a)+El{a)+2Ei{cx)E2ia)cos{a-k+->jj{t))da.{12) 
From equation (5), 



E,{a) = y/T{a) 



(13) 



for i — 1,2, where pi and p2 are the coupling efficien- 
cies functions for telescopes 1 and 2 respectively (given by 
equation (4) for non-obstructed circular pupils without 
wavefront aberrations). Therefore, 



I{t) 



Ei(t) = E, eJ'(«-'=^+"*) 



I{a) 



—pi{oL) + —P2{a) 



(6) 



for j = 1,2, with uj = 2itc/X, a. = r — ro and fci,fc2 are 
the positions of the 2 telescopes. The beam combiner (a 
fiber coupler in this case, noted X in Fig. 2) sums the 
two electric fields, with a time-dependent phase term -(/'(f) 
between the two, and the electric field at its output is 



4 ' 4 
+2^^ VPi(")P2(a) cos{a ■ k + ^P{t))dc 



(14) 



E(t) = E^{t)+E2{t)e^^^'^ 



(7) 



The intensity measured at the output of the beam com- 
biner is 



I{t) ^El+El + 2E1E2 cos{a ■ k + V'(i)), 



(8) 
(9) 



3.3. Measurements of the object's Fourier transform 
without wavefront aberrations 

Several cases are briefly discussed in the absence of wave- 
front aberrations : 

1. di^d2 = d, p(a) = po 

The coupling of the fiber is considered constant across 
the field. This corresponds to the small field approxi- 
mation of fiber interferometers, a valid approximation 
when the object is unresolved by the individual aper- 
tures. Equation (14) becomes 



T{cx.)da + TZe 



i{k) 



(15) 
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where T is the Fourier transform of X. The interferom- 
eter measures one component of the Fourier transform 
of the source's light distribution per baseline. 
2. d\= 6,1 = d 

Because the two apertures have the same diameter, 
the coupling efficiency is also the same for the two 
apertures for all values of ol. Equation (14) becomes 



m = — 



a + TZe 



where 

Ia{ot) = p{a.) X J(a). 



(16) 



(17) 




The interferometer measures here one Fourier com- 
ponent (per baseline) of an apodized image of the 
source. The apodization profile is the transmission 
profile of the fiber on the sky. 

3. di ^ d2 

The coupling efficiencies of the two fibers on the sky 
are different. Equation (14) becomes 



Fig. 3. Effect of the atmospheric turbulence on the cou- 
pling efficiency of a single-mode fiber on the sky. This 
figure shows typical instantaneous coupling efficiencies 
caused by atmospheric turbulence for two values of d/ro 
{d is the telescope diameter and ro is the Fried parame- 
ter defining the strength of atmospheric turbulence). The 
brightness scale is logarithmic from (no transmission) to 
1 (total transmission). 



j Iai{oi)da. 



(14) becomes : 

din 



m 



/(X TT f 
I{a)da + -^P2{t) I I{oc)da 



+ ^ [JIa2{cx)da]+^Tle 
where 

Iai{o!.) = pi{a) X J(a), 
for i = 1,2, and 



(21) 



Iai2(«) = \/pi{a)p2{a) X J{a). 



(19) 



(20) 



The complex visibility of the observed object, V{k), can 
be expressed as 



V{k) 



m 



(22) 



The interferometer measures a Fourier component of 
the product of the source and an apodization mask. 

The apodization mask is the square root of the product 
of the two fiber apodization profiles on the sky. 

3.4. Effect of time-variable wavefront perturbations 

Time-variable wavefront aberrations change the fiber 

coupling efficiencies of each fiber for an on-axis point 
soiirce. More seriously, the on-sky maps of the fiber 
coupling efficiency are changing in structure. In the 
following discussion, we note pi{a,t) and p2(a.t) the 
coupling efficiencies of the two fibers for the sky position 
a at time t. 



3.4.1. Observation of an object of small size 

Here, we consider an object which is not resolved by the 
individual apertures : the size of the object is much smaller 
than the diffraction limit of each telescope. The coupling 
efficiency is then constant across the object for each aper- 
ture. Because pi and p2 are not functions of a, equation 



The photometric outputs (Fig. 2), Pi{t) and P2{t), are 
measuring respectively pi{t) and p2{t). 



m)=Pi^Pi{t) ji{cx)dcx. 



(23) 



for i = 1,2, where /3i and /32 are the fraction of fiber 
flux being sent into the photometric outputs by the fiber 
splitter Yi and Y2 respectively. Qi and Q2 represent the 
flux in each fiber to the combiner (fiber splitter X). 



Qi{t) = 



X Pi{t) 



(24) 



for t = 1,2. 



I{t) = Qi{t)+Q2{t)WQi{t) X Q2{t)xTle [e^'^(*V(fe)J .(25) 

From equations (24) and (25), it appears that 1^(^)1 can 
be computed from the photometric and interferometric 
output(s) Pi{t), P2{t) and I{t). Precise measurement of 
ijj{t) is required to compute the phase of V{k). This self- 
referencing technique is used on the fiber-fed FLUOR 
beam combination unit of the IOTA interferometer. 
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3.4.2. Observation of an extended object 

If the size of the object is a significant fraction of the 
diffraction Umit of the telescopes, p is not uniform across 
the object: for a given telescope, the fiber coupling effi- 
ciency is a function of both the angular vector a and of 
time. Figure 3 shows a typical instantaneous coupling ef- 
ficiency p as a function of position on the sky a for two 
values of atmospheric turbulence. The photometric out- 
puts are: 

P,{t) = p,(a,t)J(a)da, (26) 

for « = 1,2. The interferometric output is: 



I{t) = Qi{t)+Q2{t)+^/Qi{t)Q2{t)ne 



) ial2{k) 
ial2{0) 



(27) 



where 



Iai2 (a) = /(a) Vpi(a,i)p2(a,i)- (28) 

In this case, the knowledge of I{t), Pi{t) and P2{t) is not 
sufficient to compute the visibility |T^(fe)| of the object: 
pi{oc,t) and p2{a,t) also need to be known to accurately 
interpret the observed fringe visibility. To illustrate this 
effect, the simple example of the observation of a double 
star observed with two identical telescopes is considered. 
The wavefront of each telescope is aff'ected by a random 
tip-tilt, due to atmospheric turbulence or telescope point- 
ing errors, ai and CX2 are the positions of the two stars, 
of identical brightness and both unresolved by the inter- 
ferometer. Three cases for an observation at time to are 
considered: 

(1) Telescope 1 is pointing at star 1 and telescope 2 

is pointing at star 2 {pi(ai,to) = 1, pi{a2,to) = 0, 
P2{ai,to) = and p2{,cx2,to) = 1). 

(2) Telescope 1 and 2 are both pointing halfway between 
star 1 and 2 {pi{ai,to) = pi{a2,to) = P2{cxi,to) = 
P2{oc2,to) = 0.5). 

(3) Telescope 1 and 2 are both pointing at star 1 
((pi(Q;i,io) = P2{ai,to) = 1, pi{cx2,to) = ^2(0:2,^0) = 
0). 

Pi (to) and ^2(^0) are identical in each case : it is impossi- 
ble to tell between the 3 cases from the values of the pho- 
tometric signals. In case 1 there is no coherence between 
the light in the two fibers and Iai2 = and the measured 
fringe visibility is 0, while in case 2, Xai2 = 0.5 x I and the 
measured fringe visibility corresponds to the separation of 
the 2 stars. Finally, in case 3, the observed visibility is 1, 
because the interferometer behaves just as if there was 
only one star. During real observations, the atmospheric 
seeing and pointing errors of the telescopes would induce 
variations into the measured visibilities which cannot be 
corrected by measurements of Pi and P2- 

4. Effect of atmospheric turbulence on object 
visibility measurements. 

In this section, the amplitude of the effect presented above 
is estimated for the observation of various astrophysical 



sources. In these numerical simulations, only the effect 
of atmospheric turbulence has been simulated, for vari- 
ous values of tq, the Fried parameter. This is adequate 
for ground-based interferometers without Adaptive Optics 
(such as CHARA, PTI, IOTA and NPOI). The effect of 
atmospheric turbulence on ground-based interferometers 
with adaptive optics (such as Keck, OHANA and VLTI) 
can still be obtained from this study by adopting a ro 
value which best represents the corrected wavefront: the 
value of ro that would yield the full width at half maxi- 
mum (FWHM) of the corrected long-exposure PSF if there 
was no AO system. For each telescope-!- AO system, this 
derived value of ro will decrease as a function of the bright- 
ness of the source used for wavefront sensing. 

To reproduce the effect of atmospheric turbulence, a 
simple 2-telescope fiber-fed interferometer has been sim- 
ulated. Two independent time series of 2000 wavefronts 
(d/ro = 1) were first generated, one for each telescope. 
These wavefronts are then scaled to the appropriate value 
of d/ro for the simulation. For each step, the simulator, 
written in C, first computes the on-sky fiber coupling ef- 
ficiency map for each telescope. These two maps are then 
multiplied by the intensity map of the source to yield 
the values of the photometric outputs (equation 26). The 
square root of the product of these two maps is multiplied 
by the intensity map of the source (equation 28), and the 
result is then Fourier transformed to obtain the uncor- 
rected fringe visibility, which is then corrected using the 
values of the photometric outputs (equation 27). Although 
the code can work with different telescope pupils (shape, 
diameter, central obstruction etc.), the results presented 
in this work were obtained with two circular pupils (no 
central obstruction) of identical diameter. 

4.1. One example : the observation of a double star 

Mass estimates of stars rely upon precise measurements of 
orbits of double stars. Thanks to interferometers, it is pos- 
sible to make such measurements on otherwise unresolved 

binaries (spectroscopic binaries). 

In this simulation, one of the stars is on the optical 
axis of each telescope. As can be seen on Fig. 4, when the 
second star is close to the edge of the field of view of the 
fiber (a/{X/d) is close to 1), the measured visibility is af- 
fected by the vignetting introduced by the fiber. One can 
correct for this effect because the electric field Eqi of the 
fundamental mode of the fiber is known. Atmospheric tur- 
bulence will introduce a variation (as a function of time) 
of the measured visibilities. One might think that it is pos- 
sible to correct for this effect by averaging a large number 
of measurements. Figure 5 shows that such an average 
(here, an average of 2000 visibility measurements, each 
one of those has been previously been corrected according 
to the values of the photometric outputs Pi and P2) is sig- 
nificantly different from the measurement that would be 
obtained without atmospheric turbulence. The difference 
is especially large when the visibility of the object is low 
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Fig. 4. Fringe visibilities for the observation of two stars 
of equal brightness separated by a (angular separation) 
without wavefront degradation {d/ro = 0). B is the dis- 
tance between the 2 telescopes of diameter d. A is the 
wavelength of the observation. On the right, B/d values 
for some interferometers have been indicated. 



100.0 




Fig. 5. This figure shows the absolute value of the mean 
difference between the fringe visibilities measurements 
with d/ro = 1 and the fringe visibilities measurements 
with d/ro = (Fig. 4). 



(the interferometer is resolving the double star) : the error 
can then be larger than the measured value of the visibil- 
ity. Figure 6 shows how these errors affect the measure of 
the system separation. 

For d/ro values larger than 1, the errors are almost 
independent of d/ro- This is an expected result since the 




10.0 



Fig. 6. This figure shows the relative error for measure- 
ment of the separation of the double star if the mean of the 
visibility measurements is used {d/ro = !)• The visibility 
of the double star is 90% along the dashed line. 



error is then a function of the ratio between the size of 
the object (separation between the 2 stars in this exam- 
ple) and the size of a speckle: this ratio defines the gradient 
of fiber coupling efhciency across the object, which is the 
cause of this effect. Figure 7 shows how those errors vary 
as a function of d/ro for a/{\/d) = 0.22 and B/d = 5.92. 
Errors for other values oi a/{\/d) and B/d vary in a sim- 
ilar way. It is found that the standard deviation between 
measurements is very close to the systematic error intro- 
duced by the atmospheric turbulence. This is true for all 
points except the ones close to a/{X/d) = 0.7, for which, 
although the standard deviation is relatively high, the av- 
erage visibility measurement is not significantly affected 
by atmospheric turbulence. This "equilibrium" value of 
a/{X/d) is a function of d/ro and is therefore not con- 
stant during an observation: this property cannot be used 
to hope to obtain exact visibility measurements. It is found 
that, in order to significantly decrease the visibility mea- 
surement error, d/ro needs to be smaller than 1. Existing 
Adaptive Optics systems deliver Strehl ratios of typically 
30% to 70% on bright sources in the near-infrared. This 
degree of correction, as suggested in Fig. 7, can only reduce 
the visibility measurement error by a factor of 2.5 at best. 
The adaptive optics correction is lower for fainter sources 
and in many cases. Adaptive Optics does not significantly 
reduce the visibility measurement error. Therefore, the re- 
sults presented in this study for d/ro — 1 are also valid 
for most ground-based interferometers, with or without 
adaptive optics. High-performance {Strehl > 0.9) could 
however greatly improve the quality of the measurements. 
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Fig. 7. Evolution of the systematic error in the visibil- 
ity measurement and the rms of the distribution of these 
measurements with the strength of the atmospheric tur- 
bulence ((i/ro). The Strehl, as measured from the set of 
wavefronts used for these simulations, is also plotted. The 
space between the two dotted lines represents the range 
of residual wavefront aberration delivered by current AO 
systems in the near-infrared. 



4.2. Other types of sources : Stellar diameter 
measurement, AGNs and YSOs 

Existing interferometers can measure stellar visibilities 
with an accuracy of about one percent. If the star is 
partially resolved by a telescope of the interferometer, 
it is however often impossible to reach this accuracy. 
We consider the observation of a star of diameter dstar 
with a 2-telescope interferometer. The error on the 
measurement of the fringe visibility due to the apparent 
size of the star and atmospheric turbulence have been 
computed with the same code as for the double star 
(averaging of 2000 independant photometry-corrected 
visibility measurements). Figure 8 shows how these errors 
affect the precision of the stellar diameter measurement 
(with d/ro = 1). The errors are generally smaller than 
in the case of the double star because the flux is more 
concentrated in the central region. Moreover, a very 
simple model (uniform disk) is adopted for the star and 
only one simple quantity (the diameter) is estimated. 
The errors tend to be large (more than a few percent) 
if the visibility is measured close to the first null. For 
B/d> 2, the accuracy on the stellar diameter is generally 
better than 1% if the measure is done at baselines shorter 
than the resolving baseline (for which the visibility is null). 

The effect of atmospheric turbulence on the fringe 
visibility measurements of AGNs and YSOs is hard to 
quantify precisely due to the lack of a simple and accurate 
model of light intensity distribution for those sources. 
However, the two examples studied above (double stars 
and stellar disks) provide us with some general estimate 
of the measurement errors for various sources and some 




10.0 



Fig. 8. Error in the measure of stellar diameters. B is 
the interferometer baseline, d the apertures diameter, a 
the stellar angular diameter. As for Fig. 6, the visibility 
along the dashed line is 90% and d/ro = 1. The errors 
are minimal along a line corresponding to a 40% visibility 
approximately, which explains why the 1%, 0.5% and 0.1% 
error curves are found on either side of this line. The large 
errors (more than 5%) are encountered when the fringe 
visibility is small. 

conclusions can be drawn: 



— Atmospheric turbulence seriously alters visibility mea- 
surements of sources that are not very concentrated: 
stellar diameters measurements are less corrupted than 
double star separation measurements. 

— The error on the measured physical quantity is 
larger when the visibility of the object is small. 
Measurements near or beyond the first zero of the vis- 
ibility curve should be avoided. 

— This effect still allows measurements of simple physical 
quantities (stellar diameters) but is likely to seriously 
alter more evolved measurements (small variations of 
stellar diameters, limb darkening, aperture synthesis 
imaging) . 

5. Discussion 

Imaging of extended sources with fiber-fed interferometers 
is made very delicate because of two effects. 

— The field of view of the single mode fiber is small (A/d). 

— When the source size is a significant fraction of the field 
of view of the single mode fiber, the visibility measure- 
ments are corrupted by atmospheric turbulence. 

Field of view of the fiber: To increase the field of view 

of the single mode fiber, the telescope diameter could be 
reduced, at the cost of lowering the flux in the fibers. 
This solution is of limited use because it seriously affects 
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the limiting magnitude of the interferometer. Another 
obvious solution is to increase the operating wavelength. 



Visibility measurement errors due to atmospheric turbulence: 
This effect has been studied in §4. One might think 

that because these errors can be estimated, the visibil- 
ity measurements could be corrected for, using this work. 
Unfortunately, during actual observations, d/ro does not 
stay constant and additional effects (dome seeing, mir- 
ror seeing, shape of the primary mirror etc.) prevent such 
correction: the wavefront perturbation key parameters are 
not constant in time. This effect is also a function of 
the morphology of the source, which is generally poorly 
known. However, several solutions to this problem can be 
implemented without much complexity : 

— Better statistical analysis of the measurements 
Measurements could be "weighted" according to the 
values of the photometric outputs. For example, visi- 
bility measurements are often very corrupted when at 
least one of the photometric outputs is low. 

— Reducing the size of the apertures Figures 6 and 
8 show that reducing d (the point of measure in these 
figures is then moved to the upper left, along a line 
of constant visibility, parallel to the dashed line) de- 
creases the visibility measurement error. This tech- 
nique is only efficient for bright (the apertures can be 
stopped down) sources. 

— Monitoring the PSF Real-time; imaging of the two 
PSFs on the; fiber heads can yield to an estimate of 
pi{oL,t) and f)2{a,t) (equations 26 to 28). 

The statistical analysis of the measurements can help 
with the standard deviation between measurements, but 
most of the systematic visibility measurement errors 
cannot be corrected for. Reducing the size of the aperture 
affects the limiting magnitude of the interferometer and 
degrades the precision of the measurements because of 
photon noise. Monitoring of the PSF seems to be the 
most promising technique and can be done by imaging 
the image reflected by the fiber head with a high-speed 
camera. Alternatively, efficient photon-counting detectors, 
such as APDs, could be used, since only a small number 
of "pixels" are required. This last technique is efficient 
for sources smaller than X/d, because large sources will 
tend to "blur" the PSF (the information that leads to 
estimates of pi and p2 is lost). 

One elegant solution to both effects is to use several 
fibers per telescope. By segmenting an aperture into sub- 
apertures and using one fiber per subaperture, the field of 
view of the fibers is increased and the visibility measure- 
ment errors due to atmospheric turbulence are reduced. 
This approach offers the advantage of reducing the size of 
the apertures (see above) without significant loss of light. 




Fig. 9. The telescope beam is coupled in several single- 
mode fibers. A lenslct array in the pupil plane is used to 
couple each subpupil's light in a single-mode fiber. 



6. The multifiber apertures 

6.1. The principle 

The use of one fiber per aperture limits the FOV to 
the size of a coherent domain in the focal plane, X/d. 
Increasing the FOV can be done by increasing the size 
of this coherent domain (reducing d) or by coupling N 
coherent domains in N single-mode fibers. There are 
therefore two solutions to increase the FOV. 



— Reducing the size of the pupil to be coupled in a fiber. 
To avoid losing most of the signal, the pupil is seg- 
mented into N subpupils and the light of each subpupil 
is coupled to a fiber. 

— Paving the focal plane with an array of fibers. Each 
coherent domain is coupled to a fiber. Alternatively, 
this technique could be implemented sequentially, by 
mapping the field of view with numerous pointings of 
the array. 

Paving the focal plane with fibers would increase the 

FOV but poses problems for beam recombination. We 
have therefore chosen to use the first solution. A possible 
implementation of this scheme is shown on Fig. 9: a lenslet 
array in a pupil plane is used to couple each subaperture 
to its corresponding fiber. This is very similar to the 
optical scheme used to couple subapertures to multi-mode 
fibers in curvature adaptive optics systems. 

All subapertures of the interferometer have the same 
physical size. As a result, the number of fibers on a given 
aperture is proportional to its area. 




Olivier Guyon: Wide field interferometric imaging with single-mode fibers 

(^b) Effect of a central obscuration on fiber coupling 



(C) 




M, = 



dsn. 



The field of view of the interferometer is 
A 



FOV= — 



(29) 



(30) 



for i = 1...N. At each telescope, the FOV is multiplied by 
the square root of the number of fibers used. Paving the 
focal plane with optical fibers would have resulted in the 
same increase of FOV. 

6.3. Coupling efficiency 

Figure 10 shows four possible subaperture geometries. 
We compute the cotipling efficiency for each of those 
geometries, neglecting tlic edge effects (Mj ^ 1). 

In geometry (a) and (b), the coupling efficiency is the 
product of the filling factor of the disks across the pupil 
and the coupling efficiency of a circular pupil into a fiber: 



p„ = - X 0.82 = 0.64 



Pb 



4sm(f) 



X 0.82 = 0.74. 



(31) 
(32) 



For geometries (c), and (d), we have computed the opti- 
mal coupling efficiency for square and hexagonal apertures 
using equation (1). 













Pupil surface area 






^xCoupling 




Flux coupled 


^v^^efficiency 



Central obscuration 



Fig. 10. Four possible subapertures geometries. 
6.2. Field of view 

N is the number of apertures in the array, di,i = 1...N, 
the diameter of each of those apertures, dga the diameter 
of each subaperture and Mi,i = 1...N, the number of 
fibers for each aperture. Since dsa should be the same for 
all subapertures. 



Fig. 11. The central obscuration of a pupil reduces the 
coupling efficiency. The pupil surface area is normalized to 
the unobstructed pupil surface area. The flux coupled is 
normalized to the total flux gathered by the unobstructed 
aperture. 



Pd = 0.815 



(34) 



A tight hexagonal paving of the pupil plane allows a cou- 
pling efficiency pd almost (0.5 % difFcrcncc) as good as 
a direct coupling of the entire circular pupil into a single 
fiber. For apertures with a central obstruction of diameter 
ax d, equation (2) becomes 



p{a) 



Pc = 0.793 



(33) 



/e-(-)'/o(2ra//a;2) [J^{T:dr/\f) - aJi{andr/Xf)] dr\ (35) 

Using this equation, we can compute, for each value 
of a, the combination of w, / and d which yields the 
best coupling efficiency. As can be seen on Fig. 11, the 
coupling efficiency decreases rapidly with the size of the 
central obscuration. For a central obstruction of 0.4, 
the coupling efficiency is 56%. However, if the pupil is 
subdivided into hcxagonals subapertures, each coupled to 
a fiber, the overall coupling efficiency is /94 = 0.815. 
Under diffraction-limited conditions, the division of the 
pupil in several hexagonal subpupils docs not result in a 
loss of coupling efficiency. In the case of a circular pupil 
with a central obstruction, there is a gain of coupling 
efficiency. 



6.4. Beam combining 

The high number of optical fibers in this concept is a 

strong constraint for the design of a beam combiner. Two 
possible schemes have been identified. 

— Reconstruction of the entrance pupils. The most 
intuitive approach is to optically reconstruct the in- 
coming wavefront (in fact, some approximation of it) 
for each pupil of the interferometer, as seen on Fig. 12. 
Wide-field beam recombination techniques can then be 
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Fizeau combiner pupil 




Fig. 12. One possible scheme for beam combination: the 
telescopes' pupils are optically reconstructed and fed in a 
Fizeau combiner. 



used as if each entire telescope beam had been trans- 
ported. In this scheme, the fibers (several per tele- 
scope) are only used to transport the wavefront from 
the telescopes to the beam combiner. One attractive 
approach for beam combination is a Fizeau-type beam 
combiner (shown on Fig. 12), where the beams are 
combined in a common focal plane. The pupil of the 
beam combiner can be identical to the entrance pupil. 



it can be a densified (Labeyrie 1996) version of it, or 
it can be a fixed pupil. In the first 2 cases, the beam 
combiner pupil needs to be rearranged to follow the 
geometry of the entrance pupil of the array (as seen 
by the source). 

Use of integrated optics components. A more 
promising technique is the use of a large number of 
fiber couplers to determine the complex amplitude of 
each subpupil in the array, relative to a common (for 
the whole array) phase. Phase closure techniques could 
be used between subapertures of different telescopes 
if more than 2 telescopes are us ed . Integrated optics 
components (Malbet et al 1999 , Bcrgcr et al. 1999 , 
Berger et al. 2001 ) can include a high number of 
fiber couplers in a small volume and seem to be a 
very promising technology for this beam combination 
scheme, illustrated in Fig. 13. In this case, the beam 
combiner considers the fibers as the coherent beams 
from a high number of sub-telescopes (as many as 
there are fibers). The beams of sub-telescopes belong- 
ing to the same telescope are already cophased thanks 
to adaptive optics. It would be very interesting to de- 
termine what is the optimal set of beam couplings for 
wide-field imaging with a limited number of photons. 




monopixel 
detectors 



Telescope 3 
fibers 



X-junction with 
high intersection 
angle (no coupling) 



X-junction with 
small intersection 
angle (coupler) 



waveguide 



Fig. 13. Another possible scheme for beam combination: 
multiple interferences are done between the fibers inside 
an integrated optics "optical processor" . 

6.5. Photometric calibration 

In this concept, the coupling field of view of an individual 
fiber is increased by dividing the entrance pupil into sub- 
apertures, each of those being coupled to a fiber. When 
this coupling field of view is significantly larger than the 
observed source, the effect presented in §3 and §4 is greatly 
reduced. If this coupling field of view is significantly larger 
than the atmospheric seeing, the flux in each fiber is con- 
stant in time and there is no need for photometric cali- 
bration of the interferometer. 

In practice, this second requirement is hard to meet 
because of the very large number of fibers required. 
Therefore, photometric calibration of the interferometer 
would still greatly improve the fringe visibility measure- 
ment accuracy, as it currently does for monofiber inter- 
ferometers (FLUOR instrument on the IOTA interferom- 
eter). This photometric calibration of the interferometer 
would be done by spitting the light of each fiber before the 
beam combination, yielding one photometric output per 
subaperture. Such a calibration unit could be integrated 
directly into the integrated optics processor presented in 
§6.4. It is very interesting to note that in this scheme, 
the different subapertures can be considered as different 
telescopes of an interferometric array, regardless of their 
belonging to the same main aperture. 

7. Use on a space interferometer 

7.1. (u,v) plane coverage and field of view. 

Guyon & Roddier (2001) have demonstrated that with 
only six 2-meter apertures, a full (u,v) plane coverage up 
to a 60m baseline is accessible with the rotation of an 
array. Wide field of view imaging is possible with such 
arrays if the phase information across the pupil is pre- 
served: collapsing each pupil on a single mode fiber would 
reduce the field of view. To exploit this full (u,v) plane 
coverage and reconstruct wide field images, the flux re- 
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ceived by the telescope over a large solid angle needs to 
be brought to the beam combiner. We consider such a op- 
timized array geometry and review the characteristic field 
of view quantities for a fiber interferometer using rota- 
tional aperture synthesis. In this section, we consider a 
array of N identical telescopes: do = di = ... ~ = d 
and Mo = Ml = ... = Mm = M. 

— Coupling field of view 

Using the results of §6.2, 

FOVc ^ ^ ^ Vm X 4, (36) 

dsa d 

with d, the diameter of each telescope, d^a the diameter 
of each subaperture and M, the number of fibers for 
each aperture. 

— Fourier field of view 

We first consider an array of N telescopes of infinitely 
small diameter. There are N{N — 1) + 1 (u,v) points 
measured per snapshot exposure (two symmetrical 
points per baseline plus the origin). In the case of a 
rotating array, the (u,v)-plane coverage of the final re- 
constructed image is a series of N{N—l)/2 circles (ex- 
cluding the origin). This correctly describes the (u,v) 
plane coverage if each aperture's fiux is coupled into 
a fiber (the information about the size of the pupil is 
not preserved). The maximum distance between two 
consecutive baseline lengths measured is 

d^ = e(N)— — (37) 

where e(A^) is t he efficiency of the (u,v) plane cover- 
age of the array ( |Guyon fc Roddicr 2001 ) and B is the 
longest baseline of the array. This efficiency decreases 
from 1 to 0.869 as N increases from 4 to 10. In this 
simple estimate, we will consider e{N) = 1, which is a 
valid approximation for A'^ < 10. In a (u,v)-plane cov- 
erage obtimized array, the diameter or the individual 
apertures is sufficient to fill the "gap" between consec- 
utive baselines: d = dc (Guyon & Roddier 2001). 
Reconstructing an image requires one to interpolate 
the (u,v) plane between the measured spatial frequen- 
cies: this is equivalent to a convolution in the (u,v) 
plane by a kernel of characteristic size equal to the 
distance between the known (u,v) points. In this case, 
this distance is dc- Therefore, with only one fiber per 
aperture, this imposes a limitation FOVj on the field 
of view of the reconstructed image: 



FOVi 



f 



X 



XN{N - 1) 
2B ■ 



(38) 



If no interpolation is done in the (u,v) plane, there is 
no "windowing" in the image plane, but the PSF has 
strong secondary peaks separated by FOVj: the clean 
field of view would therefore be FOVf . In the case of a 
subdivision of each pupil into M hexagonal subpupils 
(M fibers per aperture), equation (37) becomes 

B 

(39) 



dc < 



e(7V) 



N{N-l)/2 



and equation (38) becomes 
FOVf > 



2B 



(40) 



From equations (36) and (40), and noting that, for an 
array with a full (u,v) plane coverage {d — dc), with 
e{N) K. 1, equation (37) yields N'^d k. 2B, 



FOVf > FOVc. 



(41) 



Therefore, on the reconstructed image, the field of view 
will be limited by the coupling of the light into the fibers, 
and 



FOV = VM X 



A 



(42) 



Equation (42) states that the number of (u,v) points mea- 
sured is sufficient to yield a clean PSF without strong 
secondary peaks inside the field of view permitted by the 
fiber coupling. 

7.2. Spatial filtering by the fibers 

Single-mode fibers have been proposed for use on space 
interferometers because of their spatial filtering capabili- 
ties. The wavefront errors at the telescope's pupil result 
in a loss of coupling efficiency into the fiber but there is 
no loss of coherence inside the fiber: it only allows one co- 
herent mode to be coupled. Because image reconstruction 
and nulling are more tolerant to amplitude variations than 
phase variation, the constraints on the wavefront quality 
are reduced by the use of such fibers. 

In this concept, the wavefront over each subaperture is 
cleaned by this spatial filtering and the coupling efficiency 
decreases with the amount of wavefront distortion over 
this subaperture. The spatial filtering is done at a higher 
spatial frequency than if only one fiber was used per aper- 
ture. The signals in the fibers are perfectly cophased only 
if the wavefront is fiat at spatial frequencies lower than 
1/dsa- it is therefore essential to accurately control the 
quality of the wavefront at low spatial frequencies, either 
before injection into the fibers (adaptive optics) or after 
(active delay lines to cophase the fibers). Spatial filter- 
ing transforms the wavefront distortions of higher spatial 
frequencies into coupling efficiency losses, which can also 
affect the the visibility measurements. This effect can be 
corrected for by photometric calibration of the interfer- 
ometer, as discussed in §6.5. 

8. Use on a ground-based interferometer : 
OHANA 

8.1. Presentation of the OHANA project 

GHANA (Gptical Hawaiian Array for Nanoradian 
Astronomy) makes use of the existing optical and 
near-infrared telescopes on the summit of Mauna Kea 



to build a fiber interferometer (Mariotti et al. 1996 



12 



Olivier Guyon: Wide field interferometric imaging with single-mode fibers 



Mariotti et al. 199^ ). Optical fibers are preferred to clas- 
sical beam propagation for the simplicity of the beam col- 
lection at the telescope foci and the easy transport of co- 
herent light over hundreds of meters. The spatial filtering 
is also very desirable since the Strehl ratios delivered by 
the adaptive optics systems on the GHANA telescopes, in 
the near-infrared, are typically about 50% or less. 



Telescope 


D(m) 


Lat.(d) 


Long.(d) 


elev(m) 


CFHT 


3.6 


19.8252518 


155.468876 


4204.1 


Gemini 


8.1 


19.8238014 


155.469047 


4213.4 


Keck 1 


10.0 


19.8259465 


155.474719 


4159.6 


Keck 2 


10.0 


19.8265606 


155.474234 


4159.6 


Subaru 


8.2 


19.8255040 


155.476019 


4163.0 


IRTF 


3.0 


19.8262183 


155.471999 


4168.1 


UKIRT 


3.8 


19.8224315 


155.470327 


4198.5 


UH 


2.2 


19.8229911 


155.469434 


4213.6 



Table 1. Telescopes included in the simulation of the 
GHANA array, and their 3D coordinates. 



8.2. Fiber coupling Field of View 

The GHANA array, because it makes use of pre-existing 
telescopes, is made of apertures of various sizes, from the 
lOm-class telescopes (Keck I, Keck H, Subaru and Gemini) 
to the 4m-class telescopes (CFHT, IRTF and UKIRT). 
The field of view of the interferometer is the smallest FGV 
seen by the fibers. If only one fiber per telescope is imple- 
mented, the field of view is limited to 20 milli-arc-second 
at A = 1/im and 40 milli-arc-second at A = 2/im. 
A more serious problem for the imaging of extended ob- 
jects is the matching of the field of view for the different 
telescopes studied in §3.3: with only one fiber per tele- 
scope, the field of view coupled into the fibers varies by 
a factor of 2 from telescope to telescope (factor of 4 in 
area) . This is a serious problem when the observed source 
overfills the smallest field of view. Another obstacle to 
high accuracy object visibility measurements is the effect 
of residual wavefront aberrations, studied in §4. 

The use of several fibers per aperture, as presented in 
§6, would solve these two problems and allow observations 
of extended objects. Table 2 lists, for each telescope of the 
GHANA array, the approximate number of fibers required 
for an observation with a specific FGV. 



8.3. (u,v) plane coverage 

The instantaneous (u,v) plane coverage of the GHANA 
array is very sparse due to the large baselines and 
relatively small number of apertures. However, aperture 
supersynthesis greatly improves this situation: the mea- 
sured (u,v) points lie along arcs that are drawn as the 
object moves across the sky. 





0.05" 


0.1" 


0.5" 


1" 


CFHT 


1 


4 


77 


305 


Gemini 


4 


16 


377 


1505 


Keck 1 


6 


24 


588 


2351 


Keck 2 


6 


24 


588 


2351 


Subaru 


4 


16 


396 


1581 


IRTF 


1 


3 


53 


212 


UKIRT 


1 


4 


85 


340 


UH 


1 


2 


29 


114 



Table 2. Approximate number of fibers required per tele- 
scope of the GHANA array for 0.05, 0.1, 0.5 and 1 arcsec- 
ond FGV observations at 1 /im. 



The size of the telescopes of the array (up to 10 meters) 
is not negligible when compared to the distance between 
them (up to 800m). Therefore, a significant gain in (u,v) 
plane coverage can be achieved by preserving the phase 
information across the pupil. The points measured in the 
(u,v) plane are then thick arcs. Figure 14 shows the (u,v) 
coverage accessible to GHANA when observing sources at 
various declinations, when their elevation is more than 20 
degrees above the horizon. Although there is a significant 
gain in (u,v) plane coverage, the number of telescopes and 
their diameters are still too small to fill most of the fre- 
quency gaps between the arcs. The Fourier clean field of 
view, which is given by the spacing between the arcs, is 
about A/20m for GHANA, which is about twice as small 
as the field of view of the fibers of the big telescopes of the 
array. Therefore, in this case, the use of several fibers per 
aperture is mostly motivated by the need to reduce the 
turbulence-induced visibility measurement errors in the 
array and increase the coupling FGV: the size of the gaps 
in the (u,v) plane coverage is not significantly reduced by 
using several fibers per telescope. 

9. Conclusion 

Wide field imaging with fiber-fed interferometers (and 
more generally with interferometers which use spatial fil- 
tering) is made very difficult because of the small field 
of view of the fibers (or the pinhole). This work has 
also demonstrated that the spatial filtering of partially 
corrected turbulent wavefronts limits the fringe visibil- 
ity measurement accuracy. This effect is especially se- 
rious for the next generation of interferometers (Keck, 
VLTI, GHANA) which will make use of large (lOm-class) 
telescopes with adaptive optics. For such projects, this 
work has shown that photometry-calibrated fringe visibil- 
ity measurement errors will commonly be of the order 5% 
or more for extended objects (AGNs and YSGs) if only 
one fiber per telescope is used. The use of several fibers 
per aperture is an attractive solution to both improve 
the photometry-corrected visibility measurement accuracy 
and extend the field of view of these interferometers, and 
it can also allow fiber-based wide field interferometric im- 
agers to operate from space. 



Olivier Guyon: Wide field interferometric imaging with single-mode fibers 



13 




Fig. 14. GHANA maximum (u,v) coverage for objects at various declinations. 
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